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We investigate the modified gravity theories in terms of the effective dark energy models. 
We compare the cosmic expansion history and the linear growth in different models. 
We also study the evolution of linear cosmological perturbations in modified theories of 
gravity assuming the Palatini formalism. We find the stability of the supcrhorizon metric 
evolution depends on models. We also study the matter density fluctuation in the general 
gauge and show the differential equations in super and sub-horizon scales. 
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1. Introduction 

Type la supernova distance-redshift measurement shows that the present expansion 
of the universe is accelerating To explain this phenomena we need to introduce 
a homogeneous component of energy w ith a negative pressure, dubbed dark energy 
E3 or a modification of gravity EMEE] 

We can parameterize the effective dark energy component in the modified gravity 
theories The evolution equation of Hubble parameter and that of energy density 
fluctuation can be simply expressed with this parametrization. 

The Palatini formalism where the metric and the connections are treated as 
independent variables and the energy momentum tensor only depends on the metric 
leads to a different theory from what is obtained from the metric formalism. The 
Palatini formalism of f(R) gravity results in second order differential equations due 
to the algebraic relation between the curvature scalar and the trace of the energy 
momentum tensor. We will use the evolution of linear perturbations in f(R) models 
in the physical frame to check the stability of the theory. We also investigate the 
evolution of the matter density fluctuation. 
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In the next section we introduce the parametrization of modified gravities as an 
effective dark energy. We review the Palatini f(R) gravity in section 3. In section 4, 
we investigate the linear perturbation of f(R) gravity. We also derive the stability 
equation of metric fluctuations in the high curvature limit and show the stability in 
a specific model. We show the evolutions of the metric fluctuation and the density 
fluctuation in the super horizon and the subhorizon scales in section 5. We reach our 
conclusions in section 6. 



2. Modified Gravity as An Effective Dark Energy 

In general, the modified gravity can be expressed in the following form 



(1) 



where SH represents the modification to the Fricdmann equation of general relativ- 
ity and p m is the matter density. We can rewrite this equation as the effective dark 
energy term 



H 2 = 



PDE 



(2) 



where we define the effective dark energy density as 8irGpr>E = 36H. From this 
equation f2]) we can find the effective equation of state of dark energy 



1 din SH 
wde = -1 - -- 



(3) 



3 din a 

We can also express the expansion history of the universe in very simple way for 
the specific case. If SH is proportional to H (i.e. SH = aH), then we can find the 
Hubble parameter H (z) 



(4) 




Table 1. Comparison ofSH and ude m BD, DGP, metric formalism f(R) theory, and Palatini formalism f(R) models. 
F(R) = df g ^ and F is the present value of F(R). 



SH 



BD 

DGP 
Metric f(R) 

Palatini f(R) 



6 



Hi 

9 



U_ 



l ( l (Fil - /) - 3HF + 3H 2 (F - F)) 



3F 



-i + ■ 



-3H± 



\F+lHF-lF^ + ZH 2 {F -F)) -1 



l+fim 

_, 2F-2HF-4H(F -F) 

FR-f-6HF+6H 2 (F -F) 
2F-2HF-3F£^--4H(F -F) 

FR-f+3F+3HF-3F£^+6H 2 (F -F) 



In table [1] we show the expressions of the modification term to the Friedmann 
equation (SH) and the effective equation of state of dark energy (wde) for the 



I 



I 
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different modified gravity models. BD represents Brans-Dicke theory, DGP stands 
for the brane model of Dvali, Gabadadze, and Porrati^, metric f(R) means the 
f(R) theory in the metric formalism and Palatini f(R) does that in the Palatini 
formalism^. In the Palatini formalism we can relate the functions f(R) and F(R) 
to the matter energy density p m . 



Table 2. Comparisons of H(z) and evolution of the linear perturbation <5 m with the effective gravitational con- 
stant G c ff in each model. In metric and Palatini formalism the perturbation calculations are held for subhorizon 
scale (i.e. £ > H). Where Q = - ^f 2 -^ and F'(T) = ^P-. 



BD + ^^(l + z^ 8 + 2H8- AttG cSP 5 = ' 4 i — 



DGP 



2 I r H 



V(?Hk) 2+4n - (1 + 2)3 ) 5 + 2^-4^^5 = g(i+ 31i+2 ^ de ] ) 



Metric /(R) JM + Q^ 5 (1 + z) 3 <5 + 2H<5 - 47rG cf fp<5 = 2( '" 2 o 0) g 



Palatini /(ii) + ng 5 (1 + z) 3 <5 + if<5 - 47rG efr pcS = -- 



2-3Q F 

F' k 2 

4n 2F+3F' p TP 



We show the expansion rate H{z)/Hq where Hq is the present value of the 
Hubble parameter, the evolution equation for the energy density perturbation S m , 
and the effective gravitational constant G e ff for different models in table [2] The 
evolution equations of d m for the f(R) models are suitable for the subhorizon scale. 
The scale dependence on G e g is the one of the major differences compared to the 
general relativity. 



3. Palatini f(R) Gravity 

We consider a modification to the Einstein-Hilbert action assuming the Palatini 
formalism, where the metric g^ v and the torsionless connection T" are independent 
quantities and the matter action depends only on metric 



S 



(5) 



where ip are matter fields. Then the Ricci tensor is defined solely by the connection 

p pc* pa _|_ pa p/3 pa p/3 (r>\ 

lx li v ~ 1 fiu, a 1 fiQ, v T 1 a/3 1 1 av 3 W 



whereas the scalar curvature is given by i? = g^R^. Then we can find the Ricci 
tensor and the scalar curvature by using the metric relation 



3 1 



1 



Ruu = R» v + ^V^FVuF - jV^VvF - -g^u-OF 
R = R- i }-UF + i -l- 2 {dFf. 



(7) 
(8) 
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We can derive the field equation of f(R) gravity in the Palatini formalism from the 
above action (O 

F(R)R pi , - \g^f(R) = 8nGT^ , (9) 

where F(R) — df(R)/dR and the matter energy momentum tensor is given as usual 
form = — -j= * 5 ^ v ^gff"'" - • From the above equations we can get the generalized 
Einstein equation 

= 8irGT^ + (1 - F)R tw - FV V F + V M V„F + -(/ - + J^DF . 

z _r z z 

(10) 

From the above equation we can derive the modified Friedmann equations 



, SvrG ^ / „„F 3F 2 

„a SttG, , F TT F 3F 2 
-2H=— („ + ,) + --#---_, (12) 

where dot denote derivatives with respect to t. 

4. Linear Perturbation and Stability in Palatini f(R) Gravity 

The line element in the conformal Newtonian gauage is given by 



ds 2 = a 2 (r) 



1 + 2*(t, x) )dT Z + 1 - 2$(r, x) dx l dx 



(13) 



The main modifications for viable models with stable high curvature limits happen 
well during the matter dominated epoch and we can take the components of the 
energy momentum tensor as 

T° = -p(l + 5), T° = pd t q, Ti=0. (14) 

From the equation (|10|) . we can find the perturbed Einstein equation 



MO- - *Oi - J*«F - + ??WjF + J(v , V „ F > 

+ (3^ F+ 3^_3(^ JF _„ F j 4S| (15) 

where we use Sf(R) = F(R)5R. If we consider the ij-component of the perturbed 
equation, then we can find 

SF 

$-*=—, (16) 



I 



l 
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where we assume the null anisotropic stress. If we use the above equation (|16p . then 
we can express the other components of the perturbed Einstein equation (| 1 5|) 



3H 2 



2 r V 



IE1 

2~F 



If* 
2 F 5 " 



IE1 
2~F 



F 2 lH'F' 3F' H' 

1 hi* 

F 2 2 H F 2 F H 



1 H' F' 
2~H~F 



2~F 



IF 



+ 1 J* 



i? 



+ + $ + * 



3tf 2 



$" + v[r' 



\F[_ 
2 F 



(* + *) 



K 2 p 



-3: 



(17) 
(18) 



-! + — )*!>' - (3^ 



3 $ 



(6 



El 

~F 



H' N F' 



3- 



~H>~F 



H' 
'if 



3 * 



= 0. 



(19) 



where primes denote derivatives with respect to In a. To capture the metric evolu- 
tion, let us introduce two parameters as in the referenced 9 the deviation from £ 
conservation and e the deviation from the superhorizon metric evolution 



W \ H H' \aH 



Be. 



(20) 
(21) 



where we define the dimensionless quantity 



B = — — 

F H' 



(22) 
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From the above equations, we can find the expression for 9 and e, 
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K p 



H 



1/2 



lEl 
2TF 
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B 2 



El 
IF 



El 



El 

If 



B' 



El 
+ -Jp 

9 + 2 



B 

El 
IF 



El 
IT 



El 
IF 



Ell 
IF 



3 -EL ]B 

2 H 



FH 



1E1 
2~H 

H' 



* - - 



H' 



B 

(23) 

3 H' 2 
2TF 2 

1 

B 

H" 
JF 

(24) 



From equation (|23p . we can recover the conservation of Newtonian gauge when 
<f> = $ and F is a constant. 

In addition to these equations, we can find very useful equation from the struc- 
ture equation © 



5F_ 
F 



IF 
~~3~F' 



(25) 



From this equation we can find the evolution equation of matter density fluctuation 



B' H' 



H" H 



H' 



5" + 2— + —B + 2— - — + 3 \8' + —B' - 2— 



B H 



H' H 




H 



H" 

> 

' H 



,E1 

H' 



B' 2 


UE1 




F 2 ~ 


- H' 


If 


H ,,2 


H' 2 


-2E1 


H' 2 


H IF 2 


H 



EL 



(26) 



Compared with previous works we do not specify the gauge of matter density to 
solve the matter density fluctuation. 

Unstable metric fluctuations can create order unity effects that invalidate the 
backgroun d ex pansion history. We can derive the evolution equation of the deviation 
parameter 221. If we differentiate the equation (f2"4")) and consider the evolution in the 
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I 
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superhorizon scale, then we have 



B W 

2 1 B 

B H 



— 3— 



H' 



~2~tP 

4) 



B 2 



,E1 



h' 



H 
3H' 



1 



B' 2 H' , 
-2— + 2—B' 



B 2 



H 2 2 H 



B + Q' 



H" 



H 
H' 



H" 



H' 



~B 



-2— - — 



H' 



H 



H' 2 H' 
Q + 4^ + 6- 



~F(%$,Hq), 



(27) 



where we use equations (fT9|) and (|21ll and F($!,<f>) is the source function for the 
deviation e and define Q as 



H" H' 



(28) 



The above equation is different from that of the metric formalism The stability of 
e depends on the sign of the coefficient of the term proportional to e. In the metric 
formalism e is stable as long as B > 0. However, the stability is complicate and 
need to be checked for each model in the Palatini formalism. 



4.1. A particular example : f(R) — (3R n 

We demonstrate the general consideration of the previous subsection with a specific 
choice for the nonlinear Lagrangian, f(R) = (3R n , where n ^ 0,2,3. The back- 
ground is simply described by a constant effective equation of state in this model. 
The Hubble parameter scales as H 2 ~ a~ 3 /™. Then it is easy to write it with its 
derivatives in terms of In a 

H' 3 H" ( 3 \ 

fl" = _ 2^ ! ~tt = \2^J • (29) 

Here the scalar curvature is R = 3(3 — n)H 2 / '(2n). From this fact, we can also find 
the derivatives of F with respect to In a 

F' _ F" _ 3(1 - n) F" _ F'" _ f 3(1 - n) 
~F ~ F 7 _ n ' ~F ~ "F 7 ~ I n 

If we use above equations (|29| and (|30j) into (|24]) . then we find that the deviation 
from the superhorizon metric evolution is null, e = 0. 




5. Evolutions of Metric and Matter Density 
5.1. Superhorizon evolution 

We consider the metric evolution in superhorizon sized, k/(aH) -C 1. In this case, 
the anisotropy relation of the equation (|2"5|) becomes 

~ (B + A\H'q, (31) 
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where A is given by 



A = - 



B(2B^+5 



(f + 1^ + ^+4 



(32) 



From the above equations we can find the superhorizon evolution equation of $ and 
8 

&' + (# +2 ¥-b + ¥L-¥- + 4-c)&+(- + — B + — + 3-cU-0. 
I B H W H I \ B H H' 



B 



B + A 



\S" 



B' H' B'A-BA' W B 

9 i_ jd _|_ 9 

B H (B + A) 2 H B + A ' ' H' H 
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H B ~ 2 B^~ 



B' H' 2 d2 
~B + lP B 



H" H' 2 



H' 
5 H 



B 



B"A - BA" (B'A - BA') (B' + A') H' B'A - BA' ( H' 



B 



(B + A) 2 



(B + A) 2 (B + A) H (B + A) 2 \ H J B + A 



H" H" H" 2 H' 2 

_Q 4 9 I 

H H' H' 2 H 2 



W 
2 H +6 



H" W 



B) 5 = {h 



(33) 



where C is defined as 

C = 



B + A + l 



B H' H" H 

h2 — B + 2 h3 

B E H' H 



(34) 



5.2. Superhorizon evolution in a particular example 

Now we can check the evolution equations in the previous subsection in a particular 
case, f(R) ~ R n . In this case, we can simplify the following quantities 

S=2(n-1), A = -4(n-l) = -2B, = JL = 
From this, we can also simplify the evolution equations |33|) and (|33|) 



(35) 



Q — An 

y = 0, 



2n 



S' = 0. 



(36) 
(37) 



The Newtonian potential <& = constant is a solution to the equation. Also the matter 
density fluctuation has the same for as general relativity, d — constant. 



5.3. Subhorizon evolution 

For subhorizon scales where k/aH ^> l,we can find the Poisson equation from the 
equation (fl"7|) 



fc 2 ($ + 1') 



2 2 

nap , 



(38) 



I 



l 
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If we use equations (|25|) and (j38j) . then we have 

K 2 p a 2 ff 2 F'\ 5 



I FH 2 k 2 F J z 
We can find the evolution equations of <I> and S 
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(40) 

'4 



(41) 



5.4. Subhorizon evolution in a particular example 

We can use the previous relation |35|) into the evolution equations (f40|) and (J4T 



3(3 -n) 3(14n 2 + 19n- 36) 



2n 



4n 2 



$ = 



y + fc^ = o 

2n 



(42) 
(43) 



The subhorizon scale evolutions of $ and 8 show different behaviors from those of 
general relativity as expected El. 



6. Conclusions 

We investigate the different modified gravity models as an effective dark energy 
models. We can parameterize the effective equation of state of dark energy in terms 
of modified gravities. We show the evolution of energy density fluctuation compared 
to that of general relativity. 

We have analyzed the stability of metric fluctuations by checking the cosmologi- 
cal evolution of linear perturbations in Palatini f (R) gravity. We have also considered 
the matter density fluctuation in the Newtonian gauge. 

We have shown that the stability of metric fluctuations in the high redshift 
limit of high curvature is not simply expressed. We need to check each model for 
the stability. However, we have found that the deviation from the super horizon 
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metric evolution is null for a specific choice of the nonlinear Einstein- Hilbert action, 
f(R) ~ R and stability of this model is guaranteed. 

We have investigated the evolution equations of Newtonian potential and matter 
density contrast in super and sub-horizon scales. In the specific model, superhorizon 
evolutions of Newtonian potential and matter density fluctuation are same to those 
of general relativity. However, subhorizon evolutions show the different behaviors 
from the general relativity case. This will give us the method to probe the possibility 
of f(R) theory. 
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